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The results from 2.5-dimensional Particle-in-Cell simulations for the interaction of a picosecond- 
long ignition laser pulse with a plasma pellet of 50-fim diameter and 40 critical density are presented. 
The high density pellet is surrounded by an underdense corona and is isolated by a vacuum region 
from the simulation box boundary. The laser pulse is shown to filament and create density channels 
on the laser-plasma interface. The density channels increase the laser absorption efficiency and help 
generate an energetic electron distribution with a large angular spread. The combined distribution 
of the forward-going energetic electrons and the induced return electrons is marginally unstable to 
the current filament instability. The ions play an important role in neutralizing the space charges 
induced by the the temperature disparity between different electron groups. No global coalescing 
of the current filaments resulted from the instability is observed, consistent with the observed large 
angular spread of the energetic electrons. 
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The fast ignition concept, first proposed more than a decade ago^, provides a possible alternative path to achieve 
ignition in inertial confinement fusion with the potential of reducing drivers energy and increasing gains significantly. 
In this concept, the ignition is first achieved in a small region (L ~ 20/im) within a compressed pellet by a ~10 ps 
intense flux of MeV electrons. The MeV electrons are produced by the absorption of a petawatt (PW)-class laser 
at the edge of the pellet. Recent experimental results in JapanQ showed 3 orders of magnitude increase of neutron 
yield with a 0.5-PW ignition laser and consequently spurred intense new research activities in this area. One of the 
research foci is to predict the ignition laser energy needed to ignite a certain target. This energy is determined by a 
number of highly nonlinear and dynamic processes such as the laser propagation in the underdense corona, the laser 
hole-boring and electron heating near the critical surface, the transport of the energetic electrons in the overdense 
plasma, and the subsequent heating of the ignition region by these electrons. No theory exists and simulations and 
Q-r experiments are the only tools to study these processes. 

Particle-in-Cell (PIC) simulations provide a first-principle, detailed description for the plasma region of the pellet, 
which comprises the majority of the pellet except for the very dense core (electron density n e > 10 23 /cm 3 ). Since the 
PIC simulations are computationally intensive, currently it is not possible to perform full scale simulations for the 
entire pellet. Figure^shows the n» and electron temperature (T e ) information of a cryogenic target, imploded by a 25- 
kJ laser, at its peak compression [3|. The T e information is represented by the plasma parameter (number of electrons 
in a Debye cube, n\^ D ) and the ratio of the electron- ion collisionality to the electron plasma frequency, v e i/w pe . In 
the region where n e < 10 23 /cto 3 , n e \^ D » 1 which satisfies the usual definition of a plasma. Notice here n e A|j is not 
overly large either, which indicates that the discrete particle effects are important. In this region v e i/uj pe << 1, which 
shows that plasma behavior is not dominated by the collisions. All these make the PIC model the ideal simulation tool 
for this region. However, even if we were to simulate a plasma slab of 250 /zm-long and 50/im x 50^m in cross section 
and resolve the smallest relevant physical length c/uj pe (Ace = l/3c/ui pe ) for n e = 5 x 10 22 /cm 3 , the number of cells 
needed would be 1.3 x 10 12 . Assuming 4 particles/cell for each species (electrons and ions), the number of particles 
needed would be 10 13 . The total memory required to store the particle momentum and position information would be 
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560 Tera-Bytes (TB). Assuming a time step limited by Courant condition of At = 1/(3^3 V^ 1 and total simulation 
time of 10 ps, the total number of steps needed would be 6.7 x 10 5 . The total number of floating-point operations 
(FLOPs) would be 2.7x 10 20 (assuming 40 FLOPs/particle-step). This would take 10 months on a computer capable of 
10 Tera FLOPs/second (TFLOPS). This kind of simulations can be routinely performed only when 100-1000 TFLOPS 
computers with large memory are widely available. Most recent fast ignition PIC simulations used scaled-down targets 
with reduced sizes and/or in two-dimension (2D) instead of three-dimension (3D) 0,0,IE 00 Many of the important 
issues, such as hole-boring, laser absorption, and current filament at ion, have been studied. 

However, some key physics depends critically on global electron paths in the target and thus depends on the target 
size and the boundary conditions used in the simulations. For example, to absorb 30% of the energy in a 10 2Q W/cm 2 
laser for 1 ps, the thickness of a plasma slab with n e — 10 22 /cm 3 needs to be 190//m, assuming every electron in 
the slab leaves with an energy increase of 1 MeV. However, the laser cannot sweep this distance in 1 ps through 
hole-boring. Even if it could, the ion layer formed from electron depletion would pull the electrons back. The actual 
absorption is achieved through continuous replenishing of cold electrons from the bulk of the target into a narrow 
region in contact with the laser front. If the simulated target is too small and cannot provide enough needed cold 
electrons, the MeV electrons may not leave the interaction region and may be heated furthermore to achieve a different 
distribution. 

In this paper we present results from a series of 2.5D (with two space dimensions but all three components of 
particle velocities and electromagnetic fields) PIC simulations with the PIC code OSIRIS 9] in a simulation box of 
100A x 100A with a round target of 51A diameter. (A is the laser wavelength and also the basic length unit the 
simulations are scaled to.) The setup (FigEJ allows the simulation of the plasma- vacuum boundary which provides an 
important alternative pathway for the MeV electrons. The target consists of a proton-electron plasma with T e = 7.4 
keV and Tj = 1 keV. The target has a uniform density core with n e = 40n c within a diameter of 32A and a coronal 
ring with the density linearly decreasing from 40n c to in 9.5A. Here n c is the critical density which corresponds 
to 1.1 x 10 21 cm~ 3 for A = 1/im. The laser is launched from the left boundary with peak intensities of / = 10 20 and 
10 21 (l/um/A) 2 W/cm 2 and a rise time of 19 laser period, after which the laser amplitude is kept constant. The laser 
transverse profile is a Gaussian with a diffraction limited spot size of w — 7.5A. Both s- and p-polarization are used 
to infer 3D effects. A typical simulation duration is 309 laser period, which corresponds to about 1 ps for A = 1/im. 

To resolve the skin depth for n = 40n c , c/u> pe rs A/40, we use a grid size of Axw pe /c = Ayuj pe /c = 0.33 with 
12032 x 12032 grid cells. Therefore Ax « 3Xd- We reduce numerical self-heating [loj by smoothing the current. The 
time step used is Ataj pe — 0.23. We use 4 particles each for the electrons and the protons (2.4 x 10 s total particles) and 
their charges are weighted according to the initial local density. The effective finite size particle collision frequency is 
roughly v e i/u) pe m 0.035 11]. The particles are reemitted with their initial temperatures when they hit the x-boundary 
and are subject to periodic boundary condition in the y-direction. The boundary conditions for the electromagnetic 
fields are open in the laser propagation direction (x) and periodic in the transverse direction (y). A typical run 
consists of 1.6 x 10 13 particle-steps. 

The rest of the paper is devoted mainly to two issues: the laser absorption and electron heating near the laser- 
plasma interface (Sec^l and the filamentation of the electron flux inside the target (Sec^. A discussion and a summary 
are provided in Sec B 

LASER ABSORPTION AND ELECTRON HEATING 

As the laser propagates through the short underdense region here and reaches the critical surface, the most prominent 
feature observed from the simulations is the filamentation of the laser accompanied by the rippling of the target 
surface (FigEJ. The observed filament width is roughly an optical wavelength. A laser can filament due to either 
relativistic mass or ponderomotive force effects. For the parameters of the fast ignition, both mechanisms will occur on 
timescales that are very short compared to the laser pulse length. Therefore, it is not as important to study the initial 
growth phase as it is to understand the eventual nonlinear steady state. The non-relativistic theory of ponderomotive 
filamentation predicts that the mode number with the fastest growth rate increases with plasma density and the 
growth rate is higher in the direction perpendicular to the polarization than in the polarization direction [lj]. These 
conclusions are verified in Fig|2| which shows that the s-polarized cases (laser electric field E out of the simulation 
plane) display stronger laser filamentation than the p-polarized cases {E in the simulation plane). With their higher 
intensity the laser filaments dig density channels at the laser-plasma interface via the ponderomotive force (Fig|3J 
right column). These micro density channels in turn focus the laser filaments like converging waveguides which dig 
deeper channels. Therefore, the laser-plasma interface is unstable to any initial transverse modulation of the laser 
intensity, which can be viewed as filamentation in the overdense plasma 13]. 
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The density channels on the laser-plasma interface can enhance the laser absorption significantly. Heating mech- 
anisms for normal-incident lasers include the normal coupling of the laser oscillation and the plasma wave at the 
critical surface^] and the J x B heating^^. The density channels provide locations where the laser E-Held can 
be in the density gradient direction. This provides additional heating mechanisms such as the enhanced resonant 
heatingflif and Brunei heating [T^|. In our 2.5D simulations, the additional heating mechanisms are present only for 
the p-polarization cases. Therefore, we expect that the laser absoprtion is larger in the p-cases than in the s-cases. 
Figure 0^, shows the laser absorption rate for the electrons, defined as the ratio of the electron kinetic energy increase 
to the incoming laser energy within a certain time interval At, r](t) — [KE(t + At) — KE(t)]/ (P^At), for both s- and 
p-polarizations and verifies this conclusion. In addition, the p-cases also show that rj increases with time, consistent 
with the fact that as the channels becomes deeper these additional heating mechanisms become more efficient. On 
the contrary rj in the s-case remains almost flat. The absorption does not change significantly when the laser intensity 
changes from 10 20 W/cm 2 to 10 21 W/ cm 2 . However, the electron energy composition changes significantly (Fig^b). 
In the / = l0 21 W/cm 2 case, more than half of the absorbed laser energy goes to the electrons with energy greater 
than 10 MeV. These super-hot electrons may be too energetic to be stopped in the target core and be useful to the 
fast ignition. The electron spectrum in the / = 10 2Q W/cm 2 cases can be approximated by a power law, which begins 
at E ~ 0.2 MeV and falls off as £H 2 " 3 ) jiL 

The additional heating mechanisms |THTl7| allow the electrons to be accelerated directly in the laser i^-field direction. 
Thus the energetic electrons are expected to to have a large transverse momentum (P y ). This can be seen from the 
P x P y phase space plot for the electrons localized in front of the laser (FigGJi). It is also clear from the phase space 
plot that the energetic electrons do not form a beam with narrow emittance. The root-mean-square angle in the 
xy-plane for the above-1 MeV electrons is plotted in FigEt>, which also shows the large angular spread for both 
10 20 W/cm 2 and l0 21 W/cm 2 cases. The large angular spread is one of the most important characteristics of the 
energetic electron distribution in these simulations and must be taken into account in studying their subsequent 
transport and in accessing the fast ignition feasibility. 



CURRENT FILAMENT INSTABILITY 



If left alone, the laser-generated energetic electrons would deposit their energy over a large area when they reach the 
dense core region due to their large angular spread. However, the transport of the energetic electrons is not a simple 
free-streaming process but a highly nonlinear one during which their interactions with emerging magnetic fields play 
a very important role. For exam ple, we have observed that the laser-heated electrons can generate intense magnetic 
fields through the VT x Vn effect [lg on the target surface (not shown here). These magnetic fields, together with the 
surface radial electric fields, cause a fraction of the energetic electrons to move along the surface through E x £?-drift, 
which actually prevents the electrons from hitting the simulation box boundary. The motion of the energetic electrons 
that are shot forward into the interior of the target is also influenced by magnetic fields, only in this case it is the 
current filament instability that provides these fields. 

To properly study the current filament instability in the fast ignition one needs to add space charge effects due 
to temperature disparity between different groups of electrons to the standard analysis of Weibel instability |0, 2o|. 
The ~100 MA current carried by the MeV electrons in the fast ignition greatly exceeds the Alfven current limit 2 11 ] - 
which is set by the pinching force on the beam electrons from the self-magnetic field and is ~30 kA for the MeV 
electrons. A return current in the background plasma must be induced to neutralize the magnetic field of the beam 
and render the whole system nearly current neutral. Unlike the previously studied system of two identical counter- 
propagating electron clouds 20J , the forward energetic electrons have a much higher transverse temperature than the 
background plasma electrons carrying the return current. Therefore, the forward current and return current pinch 
to different degrees when filamenting and a space charge imbalance will develop. The electric field from the space 
charge is along the wave vector (or the y— direction here) and resists further filamenting. The space charge effects 
are especially important when the system is in the marginally unstable regime which is the typical case for the fast 
ignition simulations here. Furthermore, the ions will respond to the electric field and their motion becomes important 
in the marginally unstable case where the growth rate is comparable to the ion plasma frequency 8, 22]. With a 
totally different origin from that in Ref. 23], the space charge effects discussed here exist even when the system is 
stable to the two-stream instability. 

The detailed analysis of the current filament instability with the space charge effects is presented in a sep- 
arated paper ji^]. Here we illustrate the basic space charge effects on the filament instability by considering a 
plasma system with 3 species: one species of hot energetic electrons moving in the +i-direction with a drift- 
ing velocity Vdi, another species of cool return electrons moving in the — cc-direction with a drifting velocity Vd2, 
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and a cold ion species. Specifically, the electron equilibrium distribution functions are assumed to be /zo(v,x) = 
nij (2-K^/vtyiVtxi) exp(-Vy/2vt yl ) exp[-(v x - V d i) 2 /2vf xl ], where n t is the density and v tx i and v ty i are the thermal 
velocities in the x- and y-directions for each electron species, (I = 1,2). Under the relevant geometry and polarization 
v z is a negligible variable. The ion equilibrium distribution function is /30 = n^5{v x )5{v y ). Both charge and current 
densities are zero in the equilibrium, ^=1 1l n l = an d n iV<ii + "2^2 = 0. Here q\ = q% = — e are the electron charge 
and (73 the ion charge. Consequently, there is no equilibrium electric or magnetic field. 

To study the linear stability of this system against the filament mode, we assume that the mode propagation vector 
is in the y-direction, k = ky, and all perturbations have spatial and temporal dependence of exp(7i — iky). For each 
species, the perturbed distribution function fn can be found from the linearized Vlasov equation, 

/^^(E + IxB).^, (1) 

uj — kv y c ov 

where mi is the mass of the particle in the l-th species. The relevant perturbed fields are E = E x x + E y y and 
B = B z z. These can be related to the charge density p and current density j x using the Maxwell's equations. The 
linear analysis yields the following two coupled equations, 

2 2 2 2 2 tt 

^ 1 + ElS-[ 1 + 6^(6)] -% + -E l 2 %Ss 2 [i + 63(G)] = 0, (2) 

z — ' k z vf 1 lo z c z — ' w — k z c z )vf , 

1=1 tyi 1=1 v > tyl 

pE#[i + 6*(6)1 + -{i + E f P f\t^ [i + 6^(6)] - 1 - o. (3) 

1=1 W 2=1 v ' { v« ;=i 

where u)„i = Aitqfni /mi is the plasma frequency of each species, = w / (y2kvt y {) , and Z is the plasma dispersion 
function |25|. 

The analysis of the dispersion relation from Eas l2l3l reveals two regimes of instability [24] . In the relative fast 
growth regime, (m/M)kvt y 2 << 7 << kvtyi, the ions can be treated as immobile (ion mass M — ► 00 and ion plasma 
frequency lj p i = 0) and the instability threshold is 
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The third term on the right hand side of Eq0]comes from the off-diagonal term in Eas l2l3l that couples the perturba- 
tions in p and j x and therefore represents the space charge effects. This term is positive definite, indicating that the 
space charges always raise the instability threshold. If the two electron species have the same transverse temperature, 
Vt y i — vt y 2, the space charge term vanishes since the current neutrality condition implies Yl^pi^dl = 0. Therefore, 
the space charges originate from the temperature disparity between the two electron groups. In general, the space 
charge term decreases the growth rate and reduces the range of k of the instability^. 

In the slow growth regime, 7 << {m / M)kvty2, the ion response must be considered. Then the instability threshold 
turns out to be 



£ c^^M > tff + + u%. (5) 

1=1 v tyi 1=1 

Compared to Eq^ we can see in this marginal instability regime, the growth rate is so small that the ions always 
have time to react to cancel any potential space charge. Therefore, the instability threshold is the same as that when 
no space charge effect is considered |20j. 

The actual electron distribution in our PIC simulations (see Fig[5Ji) cannot be accurately approximated by two 
counter-propagating maxwellians. However, the above analysis can be extended for a general electron distribution 
in the form of f e o{Px,P y ) = Yli=i$(Px — Pdi) ex P( — P y /2Pt yl ), basically breaking the distribution into N drifting 
beamlets. The thresholds in Eas l4l5l can be modified with the sum now extended to all N beamlets. We have examined 
the electron distributions in two locations in our simulations, one in the shock region and the other in the target interior 
(shown by the small boxes in FiglfJJi), for their stability property. Specifically, we plot Upi{P"di I P?yi ~ 1)> which is 
basically Eq[5]when Vt x i and k are set to zero and w p i is neglected, for each beamlets in FigEJ; andEJi. A positive 
dot contributes to the instability and a negative one to the stability. The sum of all the dots needs to be positive for 
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the instability to occur. We found that the shock region is unstable (Fig[U) but the interior region is stable (Fig[HJi). 
Furthermore, if the space charge effect is included as in EqQJeven the shock region becomes stable. Therefore, the 
electron distributions in these fast ignition simulations are marginally unstable due to their large angular spread. The 
ions play an important role to neutralize any possible space charges. This is also supported by the observation that 
the ion density always display the same filament structure as the electrons. It is also worthwhile to point out that the 
return current beamlets (those with negative Pdi) contribute more than the forward beamlets toward the instability 
because of their relatively smaller angular spread. Any current filament stability analysis without considering the role 
of the return current would be inaccurate. 

Due to dilution of the energetic electron relative density in the interior region, the electrons there are actually 
stable to the filament instability. There the magnetic pinching force cannot overcome the thermal pressure due to 
the large transverse momentum spread. The filament structure observed there is due to the electron streaming from 
the filaments in the shock region. It is therefore no surprise that no global coalescing of these filaments into a single 
strong filament is observed here (FigEb), contrary to the previous simulations with smaller transverse box size0,[g], 
since the coalescing requires that the magnetic pinching force is larger than the thermal pressure. Even without the 
global coalescing the magnetic field from the filament instability is still strong enough (reaching ^300 MG in the 
shock region for 1/xm-laser) so that the gyro radius of MeV electrons is less than the laser wavelength. The electrons 
are confined by the self-generated magnetic field near the laser spot. 

The difference regarding the filament coalescing embodies the difference in the electron distributions between 
different fast ignition simulations. In the simulations, the electron distribution can be affected by a number of factors. 
One factor is the total heat capacity of the target, which is determined by the target size and density used in the 
simulation. Given the amount of laser energy absorbed, the fraction of energetic electrons in the distribution will 
be determined by the heat capacity. The danger of using too small a target lies in that eventually all electrons in 
the simulation acquire MeV energy, a situation totally different from the actual fast ignition scheme. Another factor 
is the recycling of the energetic electrons, which is determined by the target size and the boundary conditions. If 
the transverse box size is small and no plasma- vacuum boundary is present, the electrons can recirculate through 
the laser spot many times transversely under the periodic boundary condition before they eventually leave in the 
longitudinal direction. The plasma-vacuum interface used here combined with the surface magnetic field prevents 
such transverse recirculation. However, the longitudinal recirculation is still present but occurs later than in a smaller 
sized target. Still another factor is the true 2D nature of the simulation afforded by the large transverse size used here. 
This allows the realistic simulation of the laser filaments and density channels which are important for the heating 
mechanisms | id. fl7) which produce the large angular spread observed here. 



SUMMARY AND DISCUSSION 



In this paper the results from 2.5D PIC simulations for the interaction of a picosecond-long ignition laser pulse with 
a plasma pellet of 50-/ito diameter and 40-n c are presented. The high density pellet is surrounded by an underdense 
corona and is isolated by a vacuum region from the simulation box boundary. The laser pulse is shown to filament and 
create density channels on the laser-plasma interface. The density channels increase the laser absorption efficiency and 
help generate an energetic electron distribution with a large angular spread. The combined distribution of the forward- 
going energetic electrons and the induced return electrons is marginally unstable to the current filament instability. 
The ions play an important role in neutralizing the space charges induced by the the temperature disparity between 
different electron groups. No global coalescing of the current filaments resulted from the instability is observed, 
consistent with the observed large angular spread of the energetic electrons. 

The initial target T e used here is 7.4 keV, which is higher than the actual target bulk T e w 1 keV. This T e is chosen 
mainly to avoid numerical self-heating Recently we have implemented a second-order-spline current deposition 
scheme in OSIRIS which when combined with the current smoothing allows a grid size of Ax = 12Ad with virtually 
no self- heating. We repeat the simulations with T e = 1.1 keV and see no change in the main conclusions presented 
here. 

While the generation of the energetic electrons is studied here in considerable details, their transport from the birth 
place to n — 10 23 /cm 3 needs further study. The highly nonlinear nature of the interaction between the energetic 
electrons and emerging magnetic fields require simulations with targets of ~ lOO/im-length and up to 100n c . The 
longitudinal recycling of the energetic electrons should also be eliminated through suitable boundary conditions. Only 
then can a reliable estimate of the energy flux reaching the n = 10 23 /cm 3 surface be obtained. Fortunately, with the 
advance of computer speed and PIC code improvement this type of simulations is within reach now. 

For the very dense region of n > 10 23 /cto 3 , the collisional effects become important and the electron transport 
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and its final stopping will be affected by collisional scattering^ ■ PIC or particle-fluid hybrid type of codes need 
to include a proper collision model to correctly simulate the physics there. Eventually the simulations need to be 
compared with experiments. 
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FIG. 1: (color online) Electron density and temperature (in terms of nX% and v e i/tjJ pe ) profiles of a typical target at its peak 
compression. 
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FIG. 2: (color online) 2D PIC simulation setup. 
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FIG. 3: (color online) Laser intensity Ey (left column, in unit of (mcuj/e) 2 ) and electron charge density (right column, in unit 
of n c ) for the s- (top row) and p-polarization (bottom row) cases at t — 648/oj. Here to is the laser frequency and distance is 
in unit of c/uo. 




FIG. 4: (color online) Fraction of laser energy absorbed by the electrons for different laser intensities and polarizations (a) and 
the electron energy composition in the p-polarization case for two laser intensities at t=510 fs (b). All laser intensities and 
time are calculated assuming the laser wavelength is 1 /im. 
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FIG. 5: (color online) (a) The P x P y (in unit of mc) phase space of the electrons localized near the laser-plasma interface for 
the / = 10 20 W/cm 2 and p-polarization case (t=454 fs). (b) The rms angle (in degree) of those electrons with energy above 1 
MeV for two laser intensities (both p-polarization). All laser intensities and time are calculated assuming the laser wavelength 
is 1 fim. 




FIG. 6: (color online) (a) The electron charge density (in unit of n c ) at t — 1080c/cj. The two boxes indicate where the 
stability analysis in (c) and (d) are carried out. (b) The magnetic field B z in unit of mcu)/e, which is 107 MG for lfim- laser, 
at t = 1296c/ui. The distances in (a) and (b) are in unit of c/u>. (c) The current filament stability analysis for the shock region 
(indicated by the left box in (a)), (d) The same stability analysis for the target interior (indicated by the left box in (a)). A 
dot above the horizontal axis indicates its contribution to the instability. See the text for details. The case shown here is with 
I = 10 20 W/cm 2 and p-polarization. 



